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This paper presents the results of a numerical study of hypersonic interacting flows 
at flow conditions that include those for which experiments have been conducted in the 
Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock 
(LENS) tunnel and the ONERA R5Ch low-density wind tunnel. The computations are 
made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is 
on Mach 9.3 to 11.4 flows about flared axisymmetric configurations, both hollow cylinder 
flares and double cones. The results presented highlight the sensitivity of the calculations 
to grid resolution, provide results concerning the conditions for incipient separation, and 
provide information concerning the flow structure and surface results for the extent of 
separation, heating, pressure, and skin friction. 


Nomenclature 

C’ Chapman-Rubesin viscosity constant, 

C’ =Too/iVr*/ioo » (TWn 0 - 25 

C / friction coefficient, r w f ifi.hpcyoV^o) 

C h heat transfer coefficient, q w / {O^PocV^) 

C p pressure coefficient, (p w - Poo)/(0*5pocY^) 

d maximum double cone diameter, m 

F ratio of simulated to real molecules 

k Boltzmann constant, 1.380658 xlO” 23 JRT“ 1 

Knoo,L free-stream Knudsen number, Xoo/L 
L length from cylinder leading edge to flare, m 

Li wetted length of first cone 

Moo free-stream Mach number 

n number density, m -3 

p pressure, N/m 2 

q wall heat transfer rate, k\V/m 2 

R reattachment location, m 

Reoc u l Reynolds number, PooVoo^ / Poo 

R^oc,d Reynolds number, PocVood/ Poo 

s wetted length, m 

S separation location, m 

T temperature, K 

T* reference temperature (Eq. 6), K 

u velocity in x-direction, m/s 

Yoc free-stream velocity, m/s 

V rarefaction parameter, V = M^ 

x,y model coordinates, m 

0 l incipient separation angle, deg 

Aoc mean free path in free stream, m 

■"Senior Research Engineer, Aerothermodynamics Branch, 
Fellow AIAA. 

Copyright © 2001 by the American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United States 
under Title 17, U.S. Code. The U.S. Government has a royalty- 
free license to exercise all rights under the copyright claimed herein 
for Governmental Purposes. All other rights are reserved by the 
copyright owner. 


p viscosity, Pa*s 

p density, kg/m 3 

r viscous shear stress, N/m 2 

X hypersonic viscous interaction parameter, 

x = Mlv 

Subscripts 

w wall 

oo free stream 

ref reference 


Introduction 

H YPERSONIC separated flows produce shock/ 
shock and shock/boundary layer interactions 
that create augmented aerothermal loads and reduced 
surface control effectiveness, issues critical to hyper- 
sonic vehicle design. To enhance the understanding 
of such flows, experimental and computational studies 
have been actively promoted by the North Atlantic 
Treaty Organization (NATO) Research Technology 
Organization (RTO, formerly AGARD) for basic ax- 
isymmetric configurations and for a range of Reynolds 
numbers. The current study focuses on experimental 
conditions where separation and reattachment occur 
under laminar conditions and where the flow is as- 
sumed to be steady. Calculations are made by using 
the direct simulation Monte Carlo (DSMC) method 
of Bird 1 for both hollow cylinder-flare and double- 
cone models. The flow conditions simulated include 
those for which experiments have been conducted in 
two facilities: the ONERA R5Ch Mach 10 low-density 
wind tunnel and the Calspan-University of Buffalo Re- 
search Center (CUBRC) Large Energy National Shock 
(LENS) tunnel. 
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The hollow cylinder-flare model has a sharp leading 
edge, a cylinder 101.7 mm long, and a 30° flare. The 
author has presented extensive calculations (Refs. 2 
to 8) for the cylinder-flare model at the ONERA test 
conditions highlighting the noncontinuum and contin- 
uum aspects for the flow, the sensitivity to numerical 
simulation parameters, the agreement with measure- 
ments made in the ONERA wind tunnel, differences in 
two-dimensional (2D) and axisymmetric results, com- 
parisons with boundary layer (only cylinder or plate 
results) and Navier-Stokes solutions, and the sensitiv- 
ity of the calculations for the cylinder portion of the 
model to numerical parameters. In the present pa- 
per, results of new DSMC calculations are included for 
more energetic flows, as produced in the LENS tunnel 
for Mach 11 nitrogen flow. These calculations include 
a parametric range of conditions (free-stream density 
variations) to identify incipient separation conditions 
and three test cases for which experiments have been 
conducted 9 to measure surface heating and pressure 
distributions. 

Calculations for the double-cone models 10 are for 
flow conditions similar to the hollow cylinder-flare 
study. For the double-cone configurations investi- 
gated, the shock/shock interactions are stronger than 
those for the hollow cylinder-flare models. The first 
cone half angle is 25°, while the second cone half angle 
is 55°. The current double cone configuration produces 
strong shock interactions because the attached shock 
from the first cone interacts with the stronger bow 
shock from the second cone. Also, the outer shocks are 
modified by the separation and reattachment shocks 
where the extent of flow separation is significant for 
the combinations of model size and flow conditions ex- 
amined. Results are presented that demonstrate the 
sensitivity of the surface results to numerical parame- 
ters, Reynolds number, and flow conditions. 

Viscous Interaction and Rarefaction 
Parameters 

In the present study, the characteristic dimension se- 
lected for the hollow cylinder-flare is the length from 
the cylinder leading edge to the flare L, while that se- 
lected for the double cone is the maximum diameter 
d. For the equations in this section, L will be used 
as the characteristic length. Since moderate to low- 
Reynolds number results are normally presented as a 
function of viscous interaction and rarefaction param- 
eters, definitions of those used in the current study are 
presented. 

The free-stream Knudsen number is defined as 

Kfloo,L = Aqo jL ( 1 ) 

where L is the characteristic dimension. Also, the Kn 
can be expressed in terms of Mach and Reynolds num- 
bers with the following relation 


Kn a 0yL = 1.276^7^00 fRe^L (2) 

If one uses the boundary layer thickness S as a char- 
acteristic length, then it can be shown that 

Kriocj ~ M^yjc 1 / Re^x (3) 

which is closely related to the following rarefaction pa- 
rameter 

v = MooyJc/Re^j, (4) 

where 

C’ = pV/Poo/*oo « (Too/T*) 0 ’ 25 (5) 

and p* and p* are evaluated at a reference temperature 
given by Ref. 11 (p. 319) 

~ = 0.468 + 0.532 ^ + 0.195 A& (6) 

The 0.25 exponent for the temperature ratio in Eq. 5 
results from assuming that the viscosity is a function 
of the temperature raised to the 0.75 power. 

A V parameter is often used for correlating pressure 
and force coefficients. This parameter has also met 
with some success in correlating heat transfer mea- 
surements. 11 The V parameter is closely related to 
the Knudsen number; therefore, the larger the value 
of 17, the larger the effect of rarefaction. 

The hypersonic viscous interaction parameter \ can 
be expressed is terms of V by the following expression 

X = MlV (7) 

The induced surface pressure arising from hypersonic 
viscous interactions is often correlated as a function 
of X- Large values of \ (x > 3) correspond to 
strong interaction, while small values of \ (x < 3) 
denote a weak interaction region. It has been shown 11 
that x is an appropriate viscous correlation param- 
eter for induced pressure changes in both the weak 
and the strong interaction regions. However, as the 
flow becomes more rarefied when the leading edge is 
approached, the continuum assumptions break down, 
and the measured pressures are considerably less than 
those predicted by strong interaction theory or those 
calculated with DSMC. 12 The strong interaction 
model has a definite upstream limit that appears to 
be determined by the rarefaction parameter V rather 
than X‘ The thin Rankine-Hugoniot shock structure 
breaks down upstream of V « 0.15, where there is a 
large discrepancy between the measured and Rankine- 
Hugoniot values of shock density. (Agreement between 
experimental and DSMC results are shown to be very 
good in Fig. 12 of Ref. 12). Consequently, 0.15 < V < 
0.20 is often quoted as the boundary between the 
strong interaction and merged layer regimes. 
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As will be demonstrated later, all the experimental 
test cases considered in the present study are in the 
strong interaction regime based on the previous ex- 
pression for x when evaluated at the location of the 
deflected surfaces (flare or second cone). 

DSMC Code 

The DSMC code used in the current study is the 
general 2D/axisymmetric code of Bird. 1, 13 The molec- 
ular collisions are simulated with the variable hard 
sphere (VHS) molecular model. Energy exchange 
between kinetic and internal modes is controlled by 
the Larsen-Borgnakke statistical model. 14 For the 
present study, the simulations are performed by us- 
ing nonreacting gas models while considering energy 
exchange between translational, rotational, and vibra- 
tional modes. A constant rotational relaxation colli- 
sion number of 5 was used for the calculations. The 
vibrational collision number was 50. The reference 
conditions for the VHS model were as follows: ref- 
erence temperature = 300 K, temperature exponent 
of the viscosity coefficient = 0.75, and reference diam- 
eters for 0 2 and N 2 were 3.96 x 10“ 10 m and 4.07 x 
10“ 10 m, respectively. The model surface is assumed 
to have a specified constant temperature. Full thermal 
accommodation and diffuse reflection are assumed for 
the gas-surface interactions. 

Common to the DSMC simulations presented is the 
treatment of the computational domain griding, which 
consisted of an arbitrary number of regions. Each 
region is subdivided into cells, and the cells in se- 
lected regions are subdivided into subcells to enhance 
the spatial resolution used to select collision partners. 
In general, the cell dimensions within a region were 
nonuniform in both directions, with geometric stretch- 
ing exceeding an order of magnitude in some regions. 
Also, the macroscopic quantities are time-averaged re- 
sults extracted from the individual cells. Since the 
computational regions were not run with necessarily 
the same time step, it was essential that steady state 
conditions be established before generating the final 
time-averaged results. Steady state was assumed to 
occur when the total molecules used in the simulation, 
average molecules in each region, and surface quanti- 
ties (locations and size of the separation region, and 
heating) became essentially constant when sampled se- 
quentially over significant time intervals. 

Hollow Cylinder-Flare Results 

The ONERA hollow cylinder-flare test case con- 
sidered was formulated initially as one of the test 
problems concerning shock wave/boundary layer inter- 
actions in an AGARD Working Group 18 activity 15 
for the validation of Navier- Stokes solvers for cold 
high-speed flows where the interactions produce large 
separated regions under laminar conditions. The ini- 
tial test case generated considerable interest for code 


validation, as is evident by the AGARD activity 15 and 
several independent workshops (Ref. 16, for example). 
Interest in this problem continues with the expan- 
sion of test cases to include additional flow conditions 
and model sizes 17, 18 in the current Research Technol- 
ogy Organization (RTO, formally AGARD), working 
group 10 activities. 

A motivation for investigating these test cases with 
DSMC has been to identify the level of grid resolu- 
tion and related computational parameters that one 
must use to achieve accurate results for problems with 
complex interactions (where the grid resolution is im- 
portant in directions other than the one normal to the 
surface). Furthermore, DSMC codes provide a sim- 
ulation capability that is valid across the entire flow 
spectrum of free molecular through continuum, a desir- 
able feature for the current problems. However, there 
are practical limitations when using the DSMC ap- 
proach due to excessive computational requirements as 
one moves well into the continuum regime for multiple- 
dimensional problems. 

The current study provides numerical simulations 
for several sets of hypersonic flow conditions (Table 1), 
where some of the conditions are those for which tests 
have been conducted in the ONERA R5Ch low density 
wind tunnel and the CUBRC LENS tunnel. The LENS 
test conditions denoted as LENS- 19, LENS-11, and 
LENS-9 are conditions for which tests have been com- 
pleted by CUBRC. Results of the CUBRC tests are not 
currently available for comparison with calculations; 
however, experimental measurements for surface heat- 
ing and pressure distributions will be made available 
in Ref. 9. The remaining LENS conditions have an al- 
phabetic letter extension which denote either a pretest 
nominal condition (LENS-A) or an assumed condi- 
tion resulting from parametric variations of the free- 
stream density for the LENS-9 condition — producing 
the LENS-B through LENS-G conditions. The order 
of magnitude variation in density provides a range of 
conditions that allows one to identify the conditions 
that produce incipient separation for the current ax- 
isymmetric configuration and the nature of the viscous 
and shock/shock interactions for free-stream Reynolds 
numbers in the range of 2 500 to 25000. 

Details of the model configuration used in the ON- 
ERA tests are presented in Fig. 1. The hollow cylinder 
has a sharp leading edge with a bevel angle of 15° and 
is aligned parallel to the oncoming flow. The com- 
pression flare is inclined 30° to the cylinder and is 
terminated by a hollow cylindrical section. Additional 
information concerning model construction, materials, 
and instrumentation is given in Refs. 6 and 19. 

For the CUBRC tests, two hollow cylinder-flare con- 
figurations are used as described in Refs. 17 and 18 — 
one that mimics the ONERA model and will be re- 
ferred to as the “short” flare (current calculations use 
the dimensions shown in Fig. 1 for describing the 


3 of 27 

American Institute of Aeronautics and Astronautics Paper 2001-1027 



model outer surface for both the ONERA and the 
LENS short-flare model) and a model with a much 
longer flare. For the long-flare model, the horizon- 
tal length of the flare (referred to as the long flare) is 
118.28 mm rather than the 43.3 mm shown in Fig. 1, 
and the model terminates at the end of the flare — no 
cylindrical extension. 



Dimensions in mm 

v 



Fig. 1 ONERA hollow cylinder-flare model (x 
measured from leading edge and L = 101.7 mm). 

Calculations for R5Ch Conditions and 
Comparison with Measurements 

Figures 2 through 8 present results of the calcula- 
tions first reported in Refs. 4 and 5 that describe the 
flow-field features and surface results for the ONERA 
hollow cylinder-flare test case. The experimental value 
for free-stream Reynolds number is 18916, where the 
viscosity (3.29 x 10 -6 Pa * s) is given by the Suther- 
land expression and the characteristic dimension is the 
cylinder length L. Also presented are comparisons of 
the surface results for heating, pressure, and the extent 
of separation with the experimental measurements de- 
scribed in Refs. 6, 7, and 19. The current results are 
those obtained with the finest grid resulting from the 
grid resolution study described in Ref. 4. The pre- 
vious calculations show that the extent of separation 
is quite sensitive to the grid — -a much smaller sepa- 
ration region is obtained with a coarse grid. Data 
included in Fig. 2 provide information concerning the 
grid and simulation parameters for the finest grid so- 
lution. A four-region computational domain was used 
where each cell was subdivided into four subcells (2 
x 2). The time step in each of the four regions had 
values of 75, 75, 28, and 15 ns, respectively. 

General flow features for this test case are evident in 
Figs. 2 and 3. Figure 2 shows selected Mach contours 
and streamlines, while Fig. 3 presents the normalized 
density contours. Evident is a large separation region 
characterized by a single vortex embedded in the sub- 
sonic flow region. Calculated locations for separation 
and reattachment (denoted by S and R, respectively) 
are 76.76 mm and 134.4 mm downstream of the cylin- 
der leading edge. The shock/shock interaction occurs 


near the end of the flare where the shock layer thick- 
ness is at a minimum and the density is at a maximum, 
equal to 14.4 times the free-stream value. 

Mach 9.91 air, Re^ L =1 8,91 6 

: Fine-grid results Region Cells Subcells/Cell 
0-12 : 4-region domain 1 75x 100 2x2 

: 78,100 cells 2 100x225 2x2 

0 10 L 312,400 subcells 3 120x375 2x2 

: 1 ,900,574 molecules 4 20 x 1 80 2x2 

: At, = 75 ns 


0.08 b t, = 3.67 to 5.25 ms 



0.00 0.05 0.10 0.15 


x, m 

Fig. 2 Flow structure and simulation parameters 
for ONERA hollow cylinder flare at R5Ch condi- 
tions. 



0.00 0.05 0.10 0.15 

x, m 

Fig. 3 Density contours for ONERA hollow cylin- 
der flare at R5Ch conditions. 

Figures 4 through 6 present the calculated surface 
results for skin friction coefficient, heating rate, and 
pressure coefficient as a function of the distance from 
the cylinder leading edge (normalized by the cylin- 
der length L). Maximum values for friction, heating 
rate, and pressure occur on the flare at a location 
downstream of reattachment — very close to the end 
of the flare located at x/L = 1.426. Included in these 
figures are the results of the experimental measure- 
ments 6,7 for the extent of separation as inferred from 
oil flow measurements, heating rates extracted by us- 
ing a thin-film technique, and surface pressure inferred 
from variable reluctance differential transducers con- 
nected to the model pressure taps by tubes. 

Agreement between the calculated and the measured 
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Table 1 Free-stream and surface conditions for hollow cylinder-flare calculations. 


Condition 

too, rn/s 

Uoc, rn 

-3 

Too, K 

Poo, kg/m 3 

Poo, N/m 2 

Moo 

i,L 

Gas 

T w , K 

R5Ch 

1418.7 

8.950 x 

10 21 

51.0 

4.303 

X 

10 -4 

6.3 

9.91 

18916 

Air 

293.0 

LENS-A 

2 718.6 

1.463 x 

10 22 

194.1 

6.808 

X 

10~ 4 

39.2 

9.57 

14 920 

N-2 

297.8 

LENS-19 

2 575.6 

1.850 x 

10 22 

121.7 

8.607 

X 

10 -4 

31.1 

11.45 

25 360 

N-2 

296.1 

LENS- 11 

2 609.1 

1.089 x 

10 22 

128.9 

5.066 

X 

io- 4 

19.4 

11.28 

14490 

N-2 

297.2 

LENS-9 

2 566.4 

1.817 x 10 22 

121.1 

8.452 

X 

10~ 4 

30.4 

11.44 

24910 

n 2 

296.7 

LENS-B 

2 566.4 

1.240 x 

10 22 

121.1 

5.770 

X 

10“ 4 

20.7 

11.44 

17000 

n-2 

296.7 

LENS-C 

2 566.4 

1.021 x 

10 22 

121.1 

4.751 

X 

10~ 4 

17.1 

11.44 

14000 

N-2 

296.7 

LENS-D 

2 566.4 

8.022 x 

10 21 

121.1 

3.732 

X 

10~ 4 

13.4 

11.44 

11000 

N-2 

296.7 

LENS-E 

2 566.4 

7.293 x 

10 21 

121.1 

3.393 

X 

10“ 4 

12.2 

11.44 

10000 

N-2 

296.7 

LENS-F 

2 566.4 

5.105 x 

10 21 

121.1 

2.375 

X 

io~ 4 

8.5 

11.44 

7000 

N-2 

296.7 

LENS-G 

2 566.4 

1.817 x 

10 21 

121.1 

8.452 

X 

10~ 5 

3.0 

11.44 

2491 

N-2 

296.7 



Fig. 4 Skin friction coefficient and extent of sep- 
aration (oil-flow data from experiments). 

results are very good for the extent of separation and 
the heating rate distribution; however, there are no- 
ticeable differences for the pressure distribution. The 
separation location is the same for both sets of results. 
The calculated reattachment location occurs, however, 
somewhat forward of the experimental value, and the 
extent of the calculated separation (Ax/L) is 98% of 
the measurement. The heating-rate distributions are 
characteristic of those for laminar flows in that the 
heating shows an initial decrease at the location of 
separation (x/L = 0.76), a cusp-like behavior at the 
juncture, and a rapid increase along the flare. 

Among the surface quantities, the agreement be- 
tween the current calculations and measurements is 
the poorest for pressure. This discrepancy is particu- 
larly puzzling since the agreement for both heat trans- 
fer and the locations for separation and reattachment 
are very good. The trends of the two data sets are 
qualitatively consistent; however, the computational 
results are consistently higher than the measured val- 

5 



Fig. 5 Heat-transfer rate distributions. 

ues. As first pointed out in Ref. 4, the 42% discrepancy 
near the peak pressure location on the flare is very ob- 
vious; however, differences of this magnitude are also 
present along the hollow cylinder. In fact, this dif- 
ference is a constant. If the experimental pressure 
values are multiplied by a factor of 1.4, agreement be- 
tween the two data sets becomes very good, as shown 
in Figs. 7 and 8, where Fig. 8 presents an enlarged 
view focusing primarily on the cylinder. 

Additional results are presented in Ref. 8 that ad- 
dress the lack of agreement along the cylinder between 
calculation and measurement where the sensitivity of 
the surface and flow-field results to variations in the 
numerical parameters used in the results were dis- 
cussed earlier; in particular, to determine if the cal- 
culated surface pressure is influenced by additional re- 
finements and leading edge treatments. Surface results 
for pressure (see Ref. 8 for heating rate and flow-field 
profile sensitivity) are given in Fig. 9. The influence 
of the flare on the hollow cylinder- flare (denoted as 
H.C.F in Fig. 9) results extends slightly upstream of 
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Fig. 6 Pressure coefficient distributions. 


Fig. 8 Cylinder pressure coefficient distribu- 
tions — experimental data adjusted. 



Fig. 7 Pressure coefficients — experimental data 
adjusted. 

the x = 70 mm location, clearly evident in the surface 
pressure distribution. Also, the outflow boundary con- 
dition imposed (free stream at the end of the cylinder) 
has an influence on the surface pressure results that 
extends about 6 mm upstream of the cylinder termi- 
nation. The parametric variations for each of the five 
cylinder calculations are as follows: Case 1 was a one- 
region computational domain identical to that used 
in the H.S.F. simulation (see Fig. 2); for Case 2, the 
number of cells and subcells were increased by a fac- 
tor of three (140 x 150 cells), and the time step was 
20% of that for Case 1 (15 ns); Case 3 used a time 
step of 25 ns and the same cell resolution as H.C.F. 
and Case 1; Case 4 used a time step of 25 ns and the 
same grid resolution as the H.C.F. but had two ex- 
tra regions, with additional treatment in front of and 
downstream of the leading edge (a region upstream of 


the leading edge [-3 mm] and a region that extended 
a short distance [3 mm] downstream of the leading 
edge with five times the Ax resolution of the H.S.F 
case). Case 5 differed from Case 4 in that it included 
two additional regions to account for the influence of 
the beveled leading edge underside. Common to all 
the solutions were the four subcells/cell and a scaling 
of real-to-simulated molecules such that there were, 
on average, approximately 25 simulated molecules per 
cell. When the surface results for pressure and heating 
rate (not shown, see Ref. 8) distributions are compared 
with the H.S.F. results, there is no effect other than 
the expected results very near the leading edge -as the 
cell dimension in the x-direction (Ax) decreases near 
the leading edge, one gets an improved definition of 
the surface quantities where a local maximum occurs 
and then decreases in value as the leading edge is ap- 
proached (quantities are decreasing toward their free 
molecular values but will not achieve the free molecu- 
lar values because of upstream influence). As for the 
impact of these additional refinements on the down- 
stream flow-field quantities, no impact is evident on 
the density profiles at x/L locations of 0.3 and 0.6, as 
shown in Ref. 8. The conclusion 8 was that the DSMC 
results presented earlier for surface pressure are cor- 
rect for the cylinder and should be reasonably accurate 
for the flare, based on a constant discrepancy of 40% 
with measurements for both cylinder and flare. Also, 
the current DSMC results are in good agreement with 
those obtained by Markelov et al. 21 for the ONER A 
test case, in which a different DSMC code was used. 

Flow-field density measurements have also been per- 
formed at ONER A for the hollow cylinder-flare model 
by detecting X-ray emissions from the gas produced 
by electron beam impact. The experimental results 
have been compared 6 7 with numerical results obtained 
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Fig. 9 Effect of simulation parameters on cylinder 
pressure distributions — R5Ch conditions. 


by using the current DSMC results and two Navier 
Stokes codes. The agreement between measurements 
and calculations is somewhat mixed because neither 
DSMC nor Navier Stokes results provided consistent 
agreement with the measurements made at three lo- 
cations (x/L = 0.3, 0.6, and 0.76) along the cylinder. 
Details concerning these measurements and compar- 
isons are given in Ref. 7. 

Calculations for LENS and LENS-Like Conditions 

Four subsections are used to present the calcula- 
tions made for LENS or LENS-like flow conditions. 
The first section describes the calculations made for a 
set of nominal pretest conditions denoted as LENS- A, 
where the sensitivity of surface results and the extent 
of separation are demonstrated for grid (cell/subcell) 
resolution and length of flare (short versus long flare). 
Insight and experience gained from this portion of the 
study were then used to perform the computations for 
three CUBRC experimental test cases (LENS- 19, -11, 
and -9). Results for surface and flow-field structure are 
presented in the second subsection for each of the three 
test cases. The third subsection presents the results of 
calculations used to identify the conditions for incip- 
ient separation by arbitrarily varying the free-stream 
density of the LENS-9 test condition. The fourth sub- 
section discusses the results of the current calculations 
when presented in the form of conventional correla- 
tions. 

Tabulated results from the calculations concerning 
the extent of separation and the sensitivity of the re- 
sults to various cell/subcell combinations are included 
in Table 2, while information concerning maximum 
pressure and heating-rate values along the 30° flare 
are presented in Table 3. 


Results for LENS-A conditions pretest conditions 

This section presents results of DSMC calculations 
for flow about both the short- and long-flare hollow 
cvlinder-flare models at a nominal LENS flow condi- 
tion (Table 1), where the flow is Mach 9.57 nitrogen at 
a free-stream Reynolds number of 14 920 (characteris- 
tic length is the cylinder length, L = 101.7 mm, and 
the viscosity is given by the VHS 1 model). The flow is 
more energetic than the R5Ch conditions; however, the 
flow enthalpy is still quite low, and chemical reactions 
are neglected for this test case condition. As previ- 
ously discussed, the short flare model was assumed to 
have the same outer surface dimensions as the ONER A 
model (Fig. 1), while the long-flare model has simply 
an extended flare with the model terminated at the end 
of the flare. The long-flare model was included 17 in the 
experimental program to allow for a complete pressure 
recovery on the flare, thus making a more straightfor- 
ward comparison with theoretical models. The current 
results provide an indication of the sensitivity of the 
DSMC calculations to grid resolution, show that the 
extent of separation is much smaller than that for the 
R5Ch flow conditions, show that the results for the 
long and short flare are essentially identical within the 
domains common to the two models, and provide infor- 
mation concerning the general features of the surface 
results and flow-field structure. 
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Fig. 10 Flow structure and solution parameters 
for LENS-A condition with short flare. 

Figures 10 and 11 provide information describ- 
ing the computational domain, simulation parameters, 
flow structure (selected Mach contours and stream- 
lines), and calculated locations for separation and reat- 
tachment for the short and long flare CUBRC models, 
respectively. (The symbol F in Figs. 10 and 11 de- 
notes the ratio of real to simulated molecules.) These 
results are for the finest grid calculations resulting 
from several computations where grid refinement and 
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Table 2 Grid information, flow parameters, and separation results for hollow cylinder-flare calculations. 


Condition 

Flare 

Cells 

Subcells 

Moo 

Reoc.L 

A Uric,/, 

V 

X 

xs, mm 

mm 

Ax/L 

LENS-A 

Short 

85600 

964900 

9.57 

14 920 

0.00077 

0.064 

5.86 

91.11 

118.12 

0.265 

LENS-A 

Short 

85 600 

531 200 

9.57 

14 920 

0.00077 

0.064 

5.86 

90.70 

118.60 

0.274 

LENS-A 

Long 

142 000 

1344250 

9.57 

14 920 

0.00077 

0.064 

5.86 

91.13 

117.87 

0.263 

LENS-A 

Long 

49 625 

147500 

9.57 

14 920 

0.00077 

0.064 

5.86 

93.56 

116.09 

0.222 

LENS-A 

Long 

49 625 

49625 

9.57 

14 920 

0.00077 

0.064 

5.86 

94.79 

115.04 

0.199 

LENS-19 

Short 

88100 

979900 

11.45 

25360 

0.00055 

0.057 

7.42 

78.45 

125.68 

0.464 

LENS-19 

Short 

85 600 

85600 

11.45 

25 360 

0.00055 

0.057 

7.42 

77.61 

125.68 

0.473 

LENS-11 

Long 

142 000 

1344 250 

11.28 

14 490 

0.00094 

0.074 

9.44 

99.59 

106.35 

0.066 

LENS-9 

Long 

146 250 

1366 250 

11.44 

24910 

0.00055 

0.057 

7.40 

79.72 

125.07 

0.446 

LENS-9 

Long 

142 000 

142000 

11.44 

24 910 

0.00055 

0.057 

7.40 

77.18 

126.00 

0.480 

LENS-B 

Short 

85 600 

954 100 

11.44 

17000 

0.00081 

0.069 

9.05 

94.94 

117.84 

0.225 

LENS-C 

Short 

85 600 

954 100 

11.44 

14000 

0.00099 

0.076 

9.98 

100.86 

104.02 

0.031 

LENS-D 

Short 

85 600 

954 100 

11.44 

11000 

0.00126 

0.086 

11.21 

101.67 

102.02 

0.003 

LENS-E 

Short 

85 600 

954 100 

11.44 

10 000 

0.00138 

0.090 

11.72 

- 

- 

0.000 

LENS-F 

Short 

85 600 

436600 

11.44 

7 000 

0.00197 

0.108 

14.09 

- 

- 

0.000 

LENS-F 

Short 

85 600 

85600 

11.44 

7000 

0.00197 

0.108 

14.09 

- 

- 

0.000 

LENS-G 

Short 

60000 

60000 

11.44 

2 491 

0.00554 

0.180 

23.55 

- 

- 

0.000 

LENS-G 

Long 

96 780 

96 780 

11.44 

2 491 

0.00554 

0.180 

23.55 

- 

- 

0.000 


Table 3 Maximum pressure and heating rate values along 30° flare. 


Condition 

Flare 

x re /, mm 

P ref, N/m 2 

q ref, kW/m 2 

Pmai /P ref 

Qmax j^lref 

LENS-A 

Long 

84.15 

66.66 

22.16 

33.08 

15.05 

LENS-A 

Short 

84.15 

65.15 

21.52 

34.15 

15.61 

LENS- 19 

Short 

70.21 

63.38 

19.53 

53.27 

21.29 

LENS- 11 

Short 

92.61 

41.38 

14.63 

39.14 

16.14 

LENS-9 

Long 

71.06 

60.44 

18.93 

53.08 

20.97 

LENS-B 

Short 

83.74 

42.89 

14.92 

42.78 

17.31 

LENS-C 

Short 

90.92 

36.62 

13.65 

40.69 

15.70 

LENS-D 

Short 

90.50 

30.23 

12.31 

39.30 

14.84 

LENS-E 

Short 

90.08 

28.04 

11.95 

38.85 

14.30 

LENS-F 

Short 

90.05 

21.62 

10.43 

35.87 

13.06 

LENS-G 

Long 

87.96 

9.40 

7.12 

26.79 

9.44 


sensitivity studies were performed. Figure 12 is an ex- 
ample of the results for the long-flare model, showing 
the sensitivity of heating rate to different combina- 
tions of regions, cells, and subcells — a factor of 27 in 
sub cell resolution. The effect of grid resolution on 
heating shows the expected trend of decreased heat- 
ing with improved grid resolution outside the surface 
areas influenced by flow separation. Also, the size of 
the separation zone increases with improved grid reso- 
lution, as indicated by the tabulated results for Ax/L 
included in Fig. 12. The peak heating downstream of 
reattachment is slightly higher for the finer grid re- 
sults. 

The results presented in Fig. 12 do not demonstrate 
grid convergence; however, results for the short flare 
model, where the grid was identical to that for the 


long flare model for regions 1 through 3, showed that 
the surface results were in close agreement for two 
solutions, where one solution had the cell/subcell ar- 
rangement given in Fig. 10 (same as for the long flare 
solution for regions 1 through 3) and the one with the 
same number of cells (85 600) but only about half the 
number of subcells (531200). 

Comparisons of surface results for the short and 
long flare models are presented for heating coefficient, 
pressure coefficient, and skin friction coefficient in 
Figs. 13, 14, and 15, respectively. The results show 
that the surface quantities are essentially identical for 
surfaces common to the short and long flare models. 
For this particular LENS nominal test condition, the 
extent of separation is only 47% of that calculated 
and measured for the ONERA R5Ch test condition. 
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Fig. 11 Flow structure and solution parameters 
for LENS-A condition long flare. 
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Fig. 13 Heat-transfer coefficient results for LENS- 
A condition. 


Separation 



0.00 0.05 0.10 0.15 0.20 


x, m 

Fig. 12 Effect of grid on heating rate — long flare 
at LENS-A condition. 

With the smaller extent of separation, the calcula- 
tions suggest that the short flare is sufficiently long 
to negate any influence of the expansion at the end 
of the flare on the reattachment location. However, 
the end of flare expansion produces a thinning of the 
boundary layer and significant changes in the surface 
quantities — decreasing pressure (Fig. 14) and increas- 
ing heating (Fig. 13) and friction (Fig. 15). The same 
trends are evident for the short-flare model where the 
flow expands to the cylindrical extension at x = 0.145, 
as is clearly evident in the heating and friction results 
(also examine Fig. 16). 

Composite plots for heating rate, pressure, and skin 
friction are given in Figs. 16 and 17 for the short- and 
long-flare models, respectively. These data sets pro- 
vide information on the correlation of the three surface 



x, m 

Fig. 14 Pressure coefficient results for LENS-A 
flow condition. 

quantities and explain how T they are influenced by sep- 
aration. The general qualitative features are the same 
as previously discussed for the R5Ch test conditions. 

Figures 18 through 20 provide information describ- 
ing the flow structure for the short flare. Contours 
for nondimensional density, overall kinetic tempera- 
ture, and scalar pressure (nkT, where T is the overall 
kinetic temperature) are included. The flare-induced 
adverse pressure gradient is evident in Fig. 20 wdiere 
the isopressure lines coalesce into a separation shock 
that compresses the flow to a maximum density of 
29.6 times the free-stream value. Along the surface 
and downstream of reattachment are the locations for 
maximum density and scalar pressure, with magni- 
tudes equal to 29.6 and 57.2 times their respective 
free-stream values. When the current calculated re- 
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x, m 

Fig. 15 Skin friction coefficient results for LENS- 
A flow condition. 



Fig. 16 DSMC surface results for LENS-A condi- 
tion (short flare). 



Fig. 17 DSMC surface results for LENS-A condi- 
tion (long flare). 



Fig. 18 Density contours for LENS-A condition 
(short flare). 


suits for the LENS-A condition are compared with the 
calculated results for the R5Ch flow conditions, the 
maximum density is 3.25 times greater and occurs at 
the surface rather than in the shock layer as it does 
for the R5Ch case (Fig. 3). Also, the temperatures 
are much higher, the maximum temperature being 2.8 
times greater. For the maximum scalar pressure, the 
calculated LENS value is 6.2 times that calculated for 
the R5Ch condition. Additional details concerning the 
calculations for the LENS-A condition, particularly 
body normal profiles, are included in Ref. 8. 

The experience resulting from the calculations of the 
LENS-A condition influenced the approach that was 
used for simulating three LENS test conditions dis- 
cussed in the next section. 



Fig. 19 Overall kinetic temperature contours for 
LENS-A condition (short flare). 


10 of 27 


American Institute of Aeronautics and Astronautics Paper 2001-1027 




x, m 


Fig. 20 Scalar pressure contours for LENS- A con- 
dition (short flare). 

Calculations for LENS test cases 11, 19, and 9 

This section presents the results of DSMC calcu- 
lations for three experiments conducted at CUBRC 9 
with the LENS tunnel. The three test conditions (Ta- 
ble 1) are denoted as LENS-11, LENS-19, and LENS-9. 
Based on the \ and V parameters as previously de- 
fined, the flow for each of these test conditions is in the 
merged layer regime (V > 0.15) for at least the first 
10% of the cylinder length, while the remaining flow 
along the cylinder is in the strong interaction regime 
(X > 3). 

Even though many figures are presented describing 
the results (Figs. 21 through 58), the discussion will 
be brief since the material presented is repeated for 
each of the three test cases, and the format for presen- 
tation is similar to that used previously for presenting 
the LENS-A results. The experimental measurements 
for surface heating and pressure will be presented for 
the first time in Ref. 9. The three test conditions are 
for Mach 11 nitrogen and Reynolds numbers of 14000 
and 25000. The LENS- 11 condition is for a higher 
Mach number than the pretest LENS-A condition, but 
for approximately the same Reynolds number condi- 
tion. Consequently, the LENS- 11 viscous interaction 
and rarefaction parameters are greater than the cor- 
responding values for LENS-A (Table 2). Therefore, 
no grid sensitivity study was made for this test con- 
dition (more rarefied); instead, a grid similar to that 
used for the finest grid LENS-A calculation was used. 
As expected, results for the current calculations show 
that the flow structure for LENS-11 is different from 
that for the LENS-A condition in that a much smaller 
separation region is calculated (a Ax/L equal 0.066 
rather than 0.263), indicative of the very strong effect 
of Mach number for interacting hypersonic flows. 

The observed effect of Mach number is consistent 
with the experimental findings of Holden 22 and Del- 
ery 23 which show that the upstream influence and 
separation length decrease with Mach number for fixed 
ramp angle and Reynolds number. The LENS-11 
and LENS-A computations are for the same flare an- 
gle, essentially the same Reynolds number, and Mach 
numbers of 11.28 and 9.57, respectively. The nor- 


malized extent of separation and upstream influence 
(distance from cylinder-flare junction to the upstream 
separation location) for the two flow conditions are 
0.066/0.021 and 0.263/0.104, respectively. 

The LENS-19 and LENS-9 cases are for higher 
Reynolds number conditions than LENS- 11 and the 
flow conditions are approximately the same for these 
two test cases. The only significant difference for the 
two flow cases is with the test model configuration — 
short flare for LENS-19 and long flare for LENS-9. 
Since these two test cases were at higher Reynolds 
number conditions, two solutions were generated for 
each test case. The first solution was made with- 
out subcells, while the second solution was made w r ith 
approximately the same cell arrangement (minor tun- 
ing) and time steps, but with a significant number of 
subcells, particularly in the regions where separation 
and reattachment occurred — regions 2 and 3, respec- 
tively. The additional spatial resolution provided by 
the subcells, for the current cell resolution, has a small 
impact on the extent of separation (Table 2) and sur- 
face results (Figs. 57 and 58); in fact, the extent of 
separation is slightly smaller with sub cells, as was the 
case for the LENS-A short flare calculations. The sub- 
sequent results will show that the calculated flow-field 
and surface quantities for LENS- 19 and LENS-9 cases 
are essentially the same. Also, the calculated effect of 
short versus long flare (LENS-9 versus LENS-19) has 
no obvious impact on the flow structure or extent of 
separation (Table 2), as were the results for the LENS- 
A condition. 

Figures 21 through 53 present information for each 
of the three experimental test cases with respect to the 
following quantities: numerical simulation parameters; 
information concerning simulation history; calculated 
surface heating, pressure, and skin friction distribu- 
tions; flow-field contours for density, overall kinetic 
temperature, scalar pressure; and radial profiles for 
density, pressure, temperature, and tangential veloc- 
ity at various locations along the cylindrical portion 
of the model. These results are for the finest grid so- 
lutions (Table 2). 

The only additional information presented in this 
section that was not discussed for the LENS-A cal- 
culations is the simulation history data, Figs. 22, 33, 
and 44. These results provide a history of the simula- 
tion in terms of the total number of molecules and 
the average number of molecules per cell for each 
computational region. As noted previously, the cur- 
rent calculations are not time accurate in that the 
flow regions are not advanced with the same time 
step. Details on a regional basis are useful in identify- 
ing necessary, but not necessarily sufficient conditions 
that should be satisfied before time-averaged results 
are extracted. The current results indicate that the 
time to achieve a steady state condition increases sub- 
stantially as the extent of separation increases. This 
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behavior is demonstrated for LENS-11 and -19 (see 
Figs. 22 and 33) where the calculated extent of sep- 
aration (Ax/L) is 0.066 and 0.464, respectively, and 
the time required to approach a quasi-steady or steady 
state is much longer for the LENS-19 case. 

For the two higher Reynolds number conditions, 
nonequilibrium aspects of the flow are clearly evident 
(for example, Figs. 39 through 42 for the LENS- 19 re- 
sults) along the cylindrical portion of the model where 
temperature jump and velocity slip are present. A 
lack of translational equilibrium occurs because the 
individual translational temperature components (not 
shown) differ substantially, particularly adjacent to 
the surface. Consequently, the scalar pressure (nkT) 
of the gas adjacent to the surface is not the same as 
the surface normal pressure (p w , the normal force per 
unit area exerted by the gas on the surface). The 
scalar pressure adjacent to the surface is greater than 
the surface normal pressure. Additionally, the surface 
normal scalar pressure profiles experience substantial 
gradients adjacent to the surface along much of the 
cylinder (Figs. 40 and 51). 

Figures 54 through 56 provide a comparison of the 
calculated heating, pressure, and friction coefficients 
for the three LENS test conditions. The LENS- 19 
and LENS-9 conditions produce essentially the same 
surface results for the two test configurations — short 
and long flare, respectively. Also, the pressure over- 
expansion and recompression along the flare is more 
significant for the higher Reynolds number cases, and 
this pressure response is clearly evident in the heating 
(Fig. 54) and friction (Fig. 56) results. 

The effect of additional grid refinement by means of 
increasing the number of subcells produces very small 
changes in the extent of separation for the LENS- 19 
and -9 cases. For the calculated surface distributions, 
again, the additional grid refinement had negligible im- 
pact on the results, as is demonstrated by the heat 
transfer and pressure coefficient comparisons for the 
LENS-19 results shown in Figs. 57 and 58, respec- 
tively. Of course, if the cell resolution had been much 
coarser initially, then there would have been a very 
clear impact of additional grid refinement by either in- 
creasing the number of cells or subcell, as is somewhat 
evident in the results of Fig. 12. 
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Fig. 21 Flow structure and solution parameters 
for LENS-11 condition. 



Fig. 22 Simulation history for LENS-11 condition. 
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Region molecules/cell 




Fig. 23 Surface distributions for LENS- 11 condi- 
tion. 



Fig. 26 Scalar pressure contours for LENS-11 con- 
dition. 



Fig. 27 Scalar pressure contours along flare for 
LENS-11 condition. 



Fig. 24 Density contours for LENS-11 condition. 



Fig. 25 Overall kinetic temperature contours for 
LENS-11 condition. 
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Fig. 28 Density profiles for LENS- 11 condition. 
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Fig. 30 Overall kinetic temperature profiles for Fig. 32 Flow structure and solution parameters 
LENS-11 condition. for LENS-19 condition. 
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Region 1 time, s 

Fig. 33 Simulation history for LENS-19 condition. 



Fig. 34 Surface distributions for LENS-19 condi- 
tion. 



x, m 


Fig. 35 Density contours for LENS-19 condition. 
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Fig. 36 Overall kinetic temperature contours for 
LENS- 19 condition. 
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Fig. 37 Scalar pressure contours for LENS- 19 con- 
dition. 
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Fig. 38 Scalar pressure contours along flare for 
LENS-19 condition. 
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Fig. 39 Density profiles for LENS-19 condition. Fi S- 41 Overall kinetic temperature profiles for 

LENS- 19 condition. 




u, m/s 


Fig. 40 Scalar pressure profiles for LENS- 19 con- 
dition. 


Fig. 42 Tangential velocity profiles for LENS- 19 
condition. 
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Fig. 45 Surface distributions for LENS- 9 condi- 
tion. 


Fig. 43 Flow structure and solution parameters 
for LENS-9 condition. 



Fig. 44 Simulation history for LENS-9 condition. 



Fig. 46 Density contours for LENS-9 condition. 



Fig. 47 Overall kinetic temperature contours for 
LENS-9 condition. 
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Fig. 48 Scalar pressure contours for LENS-9 con- 
dition. 




Fig. 51 Scalar pressure profiles for LENS- 9 con- 
dition. 


Fig. 49 Scalar pressure contours along flare for 
LENS-9 condition. 




Fig. 52 Overall kinetic temperature profiles for 
LENS-9 condition. 


Fig. 50 Density profiles for LENS-9 condition. 
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Fig. 53 Tangential velocity profiles for LENS-9 Fig. 55 Pressure coefficient results for three LENS 
condition. test conditions. 




Fig. 54 Heat-transfer coefficient results for three Fig. 56 Friction coefficient results for three LENS 
LENS test conditions. test conditions. 


19 OF 27 


American Institute of Aeronautics and Astronautics Paper 2001- 1027 



0.06 r 


Cells Subcells Ax/L 



Fig. 57 Effect of grid refinement on heating rate 
coefficient for LENS- 19 conditions. 



Fig. 58 Effect of grid refinement on pressure co- 
efficient for LENS- 19 conditions. 

Conditions for incipient separation 

Calculations have been performed to identify 
the conditions for incipient separation, the largest 
Reynolds number for which the flow remains attached, 
for the current hollow cylinder-flare configuration. The 
free-stream conditions are based on the LENS-9 test 
case. The Reynolds number variation is achieved by 
varying the free-stream density, producing the condi- 
tions identified as LENS-B through LENS-G (Table 
1). Consequently, the calculations are for Mach 11.44 
nitrogen at a wall temperature of 296.7 K, and the key 
calculations are made with the short flare. Table 2 in- 
cludes information regarding the calculated locations 
and extent of separation, while Figs. 59 through 61 
present representative surface results, primarily fric- 
tion distributions. 



Fig. 59 Friction coefficient distributions — para- 
metric variations (free-stream density) to identify 
conditions for incipient separation. 



Fig. 60 Detail in separation region for various free- 
stream Reynolds numbers. 

Based on these calculations, a free-stream Reynolds 
number of 10000 is identified as the approximate con- 
dition for incipient separation. Figure 59 presents 
the skin friction coefficient distributions for several 
of the computations made for Reynolds numbers be- 
tween 2491 and 24910. The results show the expected 
trend of friction coefficient decreasing with increasing 
Reynolds number. The data also shows that the peak 
value on the flare (if the long flare had been used 
for all calculations) moves upstream with increasing 
Reynolds number. Also, the shape of the distribu- 
tions at the juncture of the cylinder and flare (x/L = 
1) is characterized by a dimple or downward point- 
ing depression that approaches zero for a Reynolds 
number of 10000. As the Reynolds number increases 


20 of 27 


American Institute of Aeronautics and Astronautics Paper 2001 - 1027 




Fig. 61 Calculated surface coefficient distributions 
for incipient separation condition. 

beyond 10000, the direction of the dimple reverses and 
points in a positive direction from a negative position 
or value. 

A more detailed view of the friction results in the 
separation region is presented in Fig. 60. Figure 61 
presents the heating, pressure and friction coefficient 
distributions for the incipient separation flow condi- 
tions. 

Note that the maximum surface values are not re- 
alized along the short flare for the lowest Reynolds 
number condition. For this flow condition, the short- 
and long-flare surface results are the same for all sur- 
face quantities except for the last few mm prior to the 
flow expansion onto the cylindrical extension. This 
rapid flow expansion produces a significant increase in 
the surface friction (not shown), as is evident in the 
other short-flare, higher Reynolds number results pre- 
sented in Fig. 59. 

Current results compared with correlations 

This section discusses how the present results com- 
pare with correlations that have evolved from analyses 
and experimental observations. 20,24,25 The first com- 
parison is with the incipient separation correlation 
based on the findings of Stollery, 25 Needham 20 and 
Holden 24 where the product of the free-stream Mach 
number and the surface deflection angle 0 X are plotted 
versus the hypersonic interaction parameter x with the 
experimental data supporting the following relation, 

Mood, = C^Z ( 8 ) 

which can also be expressed (using Eq. 7) as 

Bi = cV V (9) 

The value of C = 80 in the previous equations has 
stood the test of time in successfully correlating a sub- 
stantial body of both experimental data and numerical 


solutions. This correlation is shown as the solid line in 
Fig. 62. The correlation indicates that the domain at 
or below the solid line would be attached flow, while 
the domain above the incipient separation condition 
would experience separation. 



Fig. 62 Incipient separation as a function of the 
hypersonic viscous interaction parameter and com- 
parison with correlation (solid line) of Ref. 20. 

Also, Fig. 62 presents the current results where the 
circles denote results for attached flow and the squares 
denote separated flow. For the correlation to agree 
with current axisymmetric calculations, the value of 
the constant would have to be 100 rather than 80; 
that is, a higher incipient separation angle than the 
standard correlation. However, as has been noted by 
Stollery, 26 most of the experiments on which the cor- 
relation is based were done on low-aspect ratio flat 
plates with no side walls. A study of turbulent flows 
by Coleman and Stollery, 27 comparing 2D and axisym- 
metric tests shows an unmistakable trend in which the 
axisymmetric models have higher incipient separation 
angles. Consequently, the trend of the present find- 
ings are consistent with the trends resulting from the 
Coleman and Stollery 27 study. Also as pointed out by 
Stollery, 26 one probably would expect this trend on 
the physical grounds that the overall pressure rise for 
a given flare or flap angle is lower on the flare than on 
a 2D flap. 

Once separation occurs, the growth of the separation 
region and the state of the flow (steady and laminar 
issues), with increasing Reynolds number, is of inter- 
est. Figures 63 and 64 present the current results for 
the extent of separation as a function of Reynolds and 
Knudsen numbers, respectively. The extent of sepa- 
ration increases very rapidly with Reynolds number 
once separation is initiated. The Mach number effect 
is evident with the one data point at Mach 9.57 (long 
flare), while all remaining values are near Mach 11.4. 
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Fig. 63 Calculated extent of separation as a func- 
tion of Reoo^. 
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Fig. 64 Calculated extent of separation as a func- 
tion of free-stream Knudsen number. 

When the results are plotted as a function of Knud- 
sen number (incorporates both Reynolds and Mach 
effects), the trend is as shown in Fig. 64. Once the 
extent of separation has achieved a relatively small 
value, say 3%, the separation growth is approximately 
linear, with decreasing Knudsen number for the range 
of conditions investigated. 

A key issue with shock/shock and shock/boundary 
layer interactions is the resulting heating rate and 
pressure amplifications. Along the flare, the flow ex- 
periences a recompression shock and boundary layer 
thinning, both factors producing a peak heating that 
can be correlated with the peak pressure; i.e., 

Qmax / Qref “ A{Pmax /Pref ) (10) 

Marini 28 has presented comparisons of numerical re- 



Fig. 65 Present results compared with peak heat- 
ing correlation of Ref. 28 (solid line). 

suits with a broad range of experimental measure- 
ments for sharp leading edge ramp configurations in 
terms of Eq. 10, where a fit is obtained for A = 0.73 
and B = 0.8. The “ref” quantities are selected at a 
location upstream of the ramp or flare such that val- 
ues are unaffected by the downstream flow deflection. 
The results of Marini’s correlation are shown by the 
solid line in Fig. 65. The current results, represented 
by symbols (Fig. 65), agree well with the correlation. 
Table 3 presents detailed information concerning the 
current results — -location and magnitude of reference 
values, along with the maximum pressure and heating 
quantities along the flare. 

The current DSMC calculations for the hollow 
cylinder-flare configuration/s have provided results for 
four test conditions, one in the ONER A and three 
in the CUBRC test facilities. Comparisons with the 
ONERA data are shown to be good with the ex- 
ception of surface pressure. Additional solutions are 
presented that examined the conditions for incipient 
separation. Flow-field and surface results are pre- 
sented that demonstrate their sensitivity to numerical 
parameters and the nature of the agreement of the 
current results with existing correlations. 

Sharp Double Cone Results 

The results presented in this section are from 
Ref. 10, in which DSMC calculations were made for 
both the ONERA R5Ch and the CUBRC LENS flow 
conditions for model configurations that have been or 
will be tested. Some of these results are included 
to highlight situations where the shock/shock inter- 
actions are much stronger than those previously dis- 
cussed for the hollow cylinder-flare cases. For the 
double-cone models investigated, the first cone half 
angle is 25°, while the second cone half angle is either 
55° or 65°. The DSCM calculations for the larger half 
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angle cones have been made at R5Ch flow conditions 
for a range of model sizes and assumed variations in 
free-stream density. This range of flow and model pa- 
rameters provided both attached and separated flows, 
including large separated flow regions as is described 
in Refs. 4 and 5. For the largest Reynolds number 
flow condition (Re^d = 25 719) investigated in these 
studies, the Navier-Stokes calculations by Olejniczak 5 
indicated that flow was possibly unsteady. Recent 
calculations by Bird 30 for this problem, using a time 
accurate DSMC simulation, also show the flow to be 
unsteady with time-varying surface pressure values in 
the separated region. Since experiments have not been 
conducted for the 25/65 double cone, only a summary 
of the results for the 25/55 configuration is presented. 

Table 4 lists the free-stream conditions for the cal- 
culations. Table 5 presents information concerning the 
effect of grid resolution and flow conditions for calcu- 
lated size and location of the separation region. The 
identifiers used for the various LENS and LENS-like 
flow conditions are as follows: 1) a number extension, 
as in LENS-28, is a CUBRC run with a run number of 
28, and 2) a dual letter extension denotes a LENS-like 
condition where the first letter is the identifier used 
in Ref. 10, and the second letter “C” (cone) is used 
to avoid confusion with the identifiers used for the 
hollow cylinder-flare cases. As indicated in Table 5, 
calculations have not been made for an actual CUBRC 
test condition. The identifier code is as follows: “BC” 
denotes nominal pretest conditions for which grid sen- 
sitivity studies were made; “AC” denotes actual test 
conditions as reported in Ref. 18 (condition “A” in 
Table 1. of Ref. 18). Note that only minor differences 
exist in the flow conditions for LENS-28 and LENS- 
AC. “CC” denotes assumed conditions (adjustment of 
LENS-AC density) to achieve a specific Reynolds num- 
ber. 

Pretest (LENS-BC) and Test (LENS-AC) 
Calculations 

This section focuses on results for the LENS flow 
conditions where the maximum model diameter d is 
261.8 mm. The initial calculations were made for 
pretest nominal conditions (LENS-BC). A grid sen- 
sitivity study was conducted, and the results of this 
study, as it impacts the extent of separation (As/Li 
where “s” is surface distance measured from the cone 
vertex and L\ is the length of the first cone) and 
surface heating-rate distributions, are presented in 
Fig. 66. The qualitative characterization of the grid 
listed in Fig. 66 and Table 5 is I for intermediate, C for 
coarse, and VC for very coarse. The finest grid used in 
this exercise is described as intermediate because grid 
independence was not demonstrated — additional grid 
refinement is necessary to indicate whether the current 
intermediate grid is adequate. 

Once information became available about the test 
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Fig. 66 Effect of grid on heating- rate results for a 
pretest nominal LENS flow condition. 



Fig. 67 Flow structure and simulation parameters 
for 25°/55° double cone (LENS flow). 

conditions used in the CUBRC experiments, 18 the 
finest grid used in the pretest grid investigation was 
then used to make a calculation for the LENS-AC test 
conditions (free-stream conditions included in Fig. 68), 
and selected results are presented in Figs. 67 and 68. 
The general features of the shock layer structure are 
given in Fig. 67 where selected Mach contours are in- 
cluded, along with details of the numerical parameters 
used in the simulation. Values for the surface quan- 
tities are shown in Fig. 68. For the surface pressure 
distribution, the calculated values outside the region 
that are influenced by the shock/ boundary layer inter- 
actions are in close agreement with the inviscid cone 
values (Ref. 29) of 948 N/m 2 along the 25° cone and 
3 710 N/m 2 along the 55° cone. Details concerning 
the flow structure are presented in Ref. 10 where the 
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Table 4 Free-stream and surface conditions for sharp double cones. 


Facility 

^ oo? 

m/s 

Po o x 10 4 , 
kg/m 3 

Tloo X 10 22 , 
m -3 

Too, 

K 

Poo j 

N/m 2 

Gas 

A/oc 

T w , 

K 

CUBRC LENS-28 

2 664.0 

6.545 

1.407 

185.6 

36.05 

N-2 

9.59 

293.3 

CUBRC LENS-AC 

2 713.6 

6.808 

1.463 

194.1 

39.21 

n 2 

9.56 

297.8 

CUBRC LENS-BC 

2657.9 

7.782 

1.673 

198.9 

45.93 

n 2 

9.25 

293.0 

CUBRC LENS-CC 

2 713.6 

8.647 

1.858 

194.1 

49.80 

n 2 

9.56 

293.0 

ONERA R5Ch 

1418.7 

4.303 

0.895 

51.0 

6.30 

Air 

9.91 

293.0 


Table 5 Results of DSMC calculations for sharp double cones. 


Cone parameters 
Angles Li, mm d, mm 

Grid 

Test 

condition 


Separation 
X5, mm 

Reattach ment. 
xn, mm 

Ax/Li 

As/Li 

25°/55° 

101.6 

261.8 

I 

LENS-AC 

38340 

80.02 

99.87 

0.196 

0.265 

25 0 /55° 

101.6 

261.8 

1 

LENS-BC 

42150 

80.01 

100.10 

0.198 

0.269 

25°/55° 

101.6 

261.8 

c 

LENS-BC 

42 150 

80.60 

99.17 

0.182 

0.246 

25°/55° 

101.6 

261.8 

VC 

LENS-BC 

42150 

84.00 

97.81 

0.136 

0.186 

25°/55° 

25.8 

66.4 

I 

LENS-AC 

9 730 

21.75 

24.61 

0.111 

0.154 

25 0 /55° 

25.8 

66.4 

I 

LENS-CC 

12 360 

21.45 

24.74 

0.127 

0.179 

25°/55° 

25.8 

66.4 

I 

R5Ch 

12 360 

21.17 

25.36 

0.162 

0.229 



Fig. 68 Calculated surface results for double-cone 
model at a LENS test condition. 

calculations show that the maximum values for den- 
sity and scalar pressure are 155 and 319 times their 
respective free-stream values and that there is a maxi- 
mum overall kinetic temperature of 3 104 K. Compar- 
isons of translational and internal temperature pro- 
files show that the thermal nonequilibrium effects are 
confined primarily to the outer bow shock crossings. 
Opportunities will exist for comparing the present re- 
sults with the experimental measurements (heating- 
rate and pressure distributions) when the CUBRC 
data are released. 9 


25°/55° Cone With A Diameter of 66.45 
mm — LENS and R5Ch Conditions 

To examine the effect of different experimental test 
conditions on surface and flow-field features, a com- 
mon model (one that could be potentially tested in 
R5Ch) is used for making calculations at both ON- 
ERA R5Ch and LENS nominal test conditions. The 
model configuration is one that has been used in the 
CUBRC tests, but the maximum model diameter is 
66.45 mm rather than the 261.8-mm-diameter model 
tested. Reference 10 presents information concerning 
flow-field features and computational parameters for 
the smaller models. Figures 69 through 71 include in- 
formation concerning surface results for heating, pres- 
sure, and friction coefficients, respectively. Table 5 
provides information concerning the location and ex- 
tent of separation. 

When the calculation is made for the LENS- AC test 
conditions, the size or extent of the shock layer is 
slightly smaller, and the size of the subsonic region 
and extent of separation are noticeably reduced when 
compared with the results for the R5Ch flow condi- 
tions (see Ref. 10). Data are also included for a flow 
condition referred to as LENS-CC, where an adjust- 
ment to the density of the LENS- AC test condition is 
made to produce a larger free-stream Reynolds num- 
ber, a value equal to that of the R5Ch test condition. 
The surface heating, pressure, and friction coefficient 
results presented in Figs. 69 through 71 show to what 
extent the distributions are consistent. For the surface 
region where separation occurs, the surface coefficient 
values are much larger for the LENS test conditions 
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Fig. 69 Heat-transfer coefficient distributions for 
double-cone models. 



Fig. 70 Pressure coefficient distributions for 
double-cone models. 

than those for the R5Ch condition. Of course, the di- 
mensional heating and pressure values are much larger 
for the LENS conditions, about an order of magnitude 
difference in heating rates. The effect of increasing the 
free-stream density of the LENS- AC condition to pro- 
duce the LENS-CC conditions (same Reynolds number 
as the R5Ch case) is clearly evident in the heating and 
friction coefficient, distributions and extent of separa- 
tion (Table 5), but has a minor effect on the pressure 
coefficient distribution. 

Concluding Remarks 

Results of a computational study are presented for 
Mach 9.3 to 11.4 flow about hollow cylinder flare 
and sharp double cone models where the combina- 
tion of model configurations, size, and flow condi- 



s/L, 


Fig. 71 Friction coefficient distributions for 
double-cone models. 

tions produce a significant range of shock/shock and 
shock/boundary layer interactions. The computa- 
tions are made with the direct simulation Monte Carlo 
(DSMC) method, hence, low Reynolds number flows. 
The results presented provide insight into the nature of 
the shock interactions, their impact on surface quanti- 
ties, and the sensitivity of the results to computational 
parameters for flow conditions that can be produced 
in current ground-based facilities. 

Results of the hollow cylinder-flare calculations are 
compared with the experimental surface measure- 
ments made in the ONERA R5Ch wind tunnel (Mach 
9.91 air at a — 18916). The extent of the 

calculated separation region is very sensitive to the 
grid resolution used — a coarse grid results in a smaller 
separation region. Results for the finest grid investi- 
gated show very good agreement with the experimen- 
tal measurements for the separation and reattachment 
locations and surface heating. For surface pressure, 
the agreement between calculation and measurement 
is poor — the calculated values are uniformly high along 
both the cylinder and flare by a factor of 1.4. Addi- 
tional DSMC calculations are made for the cylinder to 
examine the impact of additional refinement of solu- 
tion parameters and leading edge treatment, and the 
results show no significant impact on the previously re- 
ported results. Based on these findings, it is believed 
the DSMC results for surface pressure are correct along 
the cylinder and should be reasonably accurate for the 
flare. 

Also, DSMC results are presented for two hollow 
cylinder-flare models, one having the same outer sur- 
face as the ONERA model and one with a longer flare, 
which has been tested for a range of conditions in 
the CUBRC LENS tunnel. Calculations are made for 
three of the LENS test conditions for nominal Mach 
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11.4 nitrogen flow at Re of 14490 to 25 360). Addi- 
tional calculations were made at lower Reynolds num- 
ber values (assumed density variations) to calculate 
the conditions for incipient separation. Information 
concerning the effect of grid resolution is presented, 
along with detailed data concerning surface results and 
flow structure. Comparisons of the calculated data 
with existing correlations based on experimental mea- 
surements are consistent with respect to the following: 
1) qualitative trends as to the effect of Mach number 
and Reynolds number on the extent of separation and 
upstream influence of a compression flare (ramp); 2) 
qualitative and quantitative agreement of peak heat- 
ing and pressure along flare; and 3) the conditions for 
incipient separation. 

For the 25°/55° double-cone configuration models, 
the extent of separation, as a function of flow condi- 
tions, is demonstrated for both LENS and R5Ch flow 
conditions. Results of a grid sensitivity investigation 
are discussed, and surface results are presented for flow 
conditions (Mach 9.56 nitrogen at a Re^^ = 38340) 
that are approximately the same as those of a LENS- 
28 test that has been completed. 

Opportunities will exist for comparing some of the 
current results for both the hollow cylinder-flare and 
double cone configurations with experimental mea- 
surements for surface heating and pressure distribu- 
tions. 
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